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Abstract 

Starting from the Hofstadter butterfly, we define lattice versions of Landau 
levels as well as a continuum limit which ensures that they scale to contin- 



^ . uum Landau levels. By including a next-neighbor repulsive interaction and 

^^ . projecting onto the lowest lattice Landau level, we show that incompressible 

Q . ground states exist at filling fractions, v = 1/3, 2/5 and 3/7. Already for 

On ■ values of lo/a ~ 2 where Iq (a) is the magnetic length (lattice constant), the 



lattice version of the u = 1/3 state reproduces with nearly perfect accuracy 



1-^ ■ the the continuum Laughlin state. The numerical data strongly suggests that 

a: 

O ■ at odd filling fractions of the lowest lattice Landau level, the lattice constant 



is an irrelevant length scale. We find a new relation between the hierarchy of 
incompressible states and the self-similar structure of the Hofstadter butterfly. 



There has been a considerable amount of numerical work on the fractional quantum Hall 
effect [0 (FQHE) for continuum models on planar or spherical geometries. The question we 
wish to address in this paper is: can one observe the FQHE on a lattice and if so, how does 
it compare to continuum results? Since on the lattice rotational invariance is broken, the 
pseudopotential |0] approach to the FQHE is not applicable and there is a priori no reason to 
expect that ground states at certain filling fractions of Hofstadter [^j bands are characterized 
by a specific power law decay of the density- density correlation (7(|r|) for small values of |f^. 
On the other hand there is an efficient way to construct incompressible quantum liquids 
(IQL) states which implicitly assumes the existence of Laughlin states 0| on a lattice: the 
parton approach @J^. The Hamiltonian which has parton wave functions as mean- field 
solutions can be derived from first principles only on a lattice 0. Hence, this approach 
requires that properties of the Laughlin states survive discretization of the plane. We show 
that even for small lattice sizes there exist incompressible states at odd denominator filling 
fractions of the lowest Hofstadter level if it has the Landau level degeneracy. These states 
reproduce with remarkable accuracy the relevant properties of the Laughlin states. 

Our starting point are fully spin-polarized electrons on a square lattice in a magnetic 
field: 

H, = -tj:cicje.pi^J!Adr\. (1) 

Here, the sum runs over next-neighbors, cl creates a spinless electron on site i, n^ = d-c^^ 
and in the Landau gauge A{x) = —B{x2,0,0). The square lattice lies in the el, 62 plane 
and $0 is the flux quantum. In order to put the system on a torus, we impose the boundary 
conditions: cl =. = exp (^Xl(O) c^- Here, L is a vector with norm equal to the linear size 
of the square lattice oriented along the Ci or 62 direction, and A{x + L) = A{x) + vXl(^)- 
To obtain a single-valued wave function the total flux traversing the lattice has to be an 
integer multiple of the flux quantum |^ . The spectrum of the above Hamiltonian yields the 
Hofstadter butterfly ||^. To make contact with the corresponding continuum Hamiltonian, 
we define the field operators cl = a J dfS{r — -Rj)\l'^(r). Expanding the phase factors in 
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equation (0) and the delta-function up to second order in the lattice constant, one obtains 
up to a constant term: Hq -^ ^-^ J dr^'< (r) (^P + lA{r)^ ^(r). Setting ta'^/h'^ = 1/2M 
where M denotes a constant mass, yields the desired continuum limit. The continuum 
Hamiltonian has an energy spectrum described by Landau levels: En = hujc{n + 1/2) where 
ujc = eB/Mc. Each Landau level has a macroscopic degeneracy given by iV$ = L^B/^o and 
the magnetic length is given by /q = ( ^^ ) 

To achieve such a situation on the lattice, we set 

p Ba? ^ r ^r ^T . . /^\ 

- = = for a A^ X A^ lattice. (2) 

q $0 niN 

Since the number of flux quanta traversing the lattice, A^$ = N/m, has to be an integer, 
A^ has to be divisible by m. Before arguing why this choice of p/q is suitable, let us note 
that in order to obtain a well defined continuum limit for a fixed value of m, we have to 
set B or equivalently Iq to a constant. This yields: a ~ 1/VN as well as L ~ ^/N. Hence, 
our continuum limit corresponds to a combined continuum (a -^ 0) and thermodynamic 
{L —>■ oo) limit. Since we want to keep the filling fraction u = Np/N<j, constant during the 
scaling to the combined continuum and thermodynamic limit, Np/N"^ ~ 1/A^. Here, Np 
denotes the number of particles on the lattice. 

For the above choice of p/q and in the Landau gauge, the boundary conditions are: 
^l+Ndi ~ ^1 ^^"^ ^'i+Nd2 ~ exp (— 27r2x/m) c^. Here, i = a{x,n). We define the partial 
fourier transform 7^ = -^ J2x ^~^^^(^a{x n)i ^ ~ 2vrnfc/A^ where n^ = 1 ... A^ to write the 
Hamiltonian (|l]) as: 



ff„ = E E (-2* ™<|^ + *) -^r - 1 {lY'Wr' + h.c) 

^Un+N) ^ ^t>)^ ^j^gj.g ^^ ^ 27r/m. (3) 



The boundary condition in the n index of the 7 operators mixes m fc-sectors. Hence, the 
resulting Hamiltonian is block diagonal and contains Nq, = N/m blocks, each of of dimension 
rriN X mN. We label the blocks with the quantum number ko, and denote the mN x mN 
matrices by A{ko). For this special choice of p/q (H), one may show that the matrices 



A{ko) are equivalent so that the spectrum of A{ko) is independent of ko. Therefore, each 
energy level of the Hamiltonian (|^) is at least A^$-fold degenerate. Numerically, one may 
see that the eigenvalues of one A{p) matrix are up to accidental degeneracies at zero energy, 
non-degenerate. The spectrum of the Hamiltonian (^) is thus described by Nm, A^$-fold 
degenerate energy levels. We identify those energy levels to lattice Landau levels (ILL). To 
proceed further, we have to diagonalize numerically the matrices A{ko) so as to obtain the 
single particle wave functions spanning a given ILL. For more general values of p/q = —N, 
the Hamiltonian again splits into N/m block diagonal matrices of dimension mN x mN. 
For n > 1 these matrices are not equivalent, and the spectrum is better described as mN, 
A^$/n-fold degenerate bands. One may see numerically that as the lattice size is increased, 
n low- lying adjacent bands merge very rapidly together to form a Landau level 0. 

Having defined ILL, one may now compare their properties with those of the continuum 
Landau levels. Consider the magnetic length which is given in the continuum by: /q = 
/^B \ ^ Pqj, q^j, choice oip/q the former equation is equivalent to u/Iq = J27c/(mN). In 
order to obtain an estimate of the magnetic length on the lattice, one notes that for a filled 
first Landau level (i.e. z/ = 1) 0: 

/N \~^ - - ( r^\ 

\-^) E<^W'^(^ + ^)>^l-expl-^l. (4) 

In order to obtain the magnetic length on a finite lattice, we have fitted our lattice-data 
for the density-density correlations (Ihs of the above equation) to the continuum form (rhs 
of the above equation). We consider two directions: r = r(l, 1) and r = r(l,0). For both 
directions, and already for small lattice sizes (i.e. N = 12) the data fits the above form 
extremely well for all considered values of m [^ . Figure 1 plots the lattice magnetic- length 
for both considered directions (symbols) as well as its value predicted by the continuum 
equations (solid line). As may be seen, the agreement, already for small lattice sizes is 
very good. It is to be noted that the rotational invariance present in the continuum is 
violated by the presence of the square lattice. However, this symmetry is restored since 
as may be seen the discrepancy between the two lattice values of /q vanishes rapidly. To 



understand this striking result we have considered the continuum eigenfunctions of Hq and 
the corrections due to finite lattice constant and volume. We find that the decay of the 
eigenfunctions remains gaussian in any order of perturbation theory. Similarly, one may 
check that the cyclotron frequency scales rapidly to its continuum value. Hence, the above 
defined ILL, reproduce very accurately the physics of the continuum Landau levels already 
on small lattice sizes. 

We show that IQL states exist on the lattice by considering a next-neighbor density- 
density interaction: 



Hi = — Y. ■ ^t^j -^ ^- j dr : p(f) A^-p(f) : (5) 

Here, : O : denotes the normal ordering of the operator O. A similar Hamiltonian was 



considered in reference |jTT| in the context of anyon superconductivity. For this continuum 
interaction (rhs of the above equation), it has been shown that Laughlin wave function 
is exact [|12|. Therefore, we should ultimately reproduce the physics of the Lauglin wave 
function. 

In order to carry out the calculations, we project the Hamiltonian onto the lowest ILL. 
This projection is justified in the limit B -^ oo. The resulting Hamiltonian is an effective 
one dimensional Hamiltonian which one may solve numerically by using a standard imple- 
mentation of the Lanczos method. The projected Hamiltonian has the same symmetries 
(conservation of the one-dimensional total momentum) as the Hamiltonian considered by 
Yoshioka et al. |T^. The scaling of the hopping matrix element, t, and interaction strength, 
V, follow from the continuum limit: V = VqN"^ and t = IqN where Vq and to ^^re constants 
independent of the lattice size. The former follows from the fact that a ~ 1/yN . 

For the numerical calculations we chose m = 1 so that A^ = Nq>. This choice gives 
the smallest magnetic length, and hence lattice effects should be maximized. The following 
refers to this special case. Our numerical analysis is based on the calculation of the quantity 
AEq{Np/N) = Eo{Np/N) - Npef^ where Eo{Np/N) is the ground state energy of the 
projected Hamiltonian for an A^-site chain with Np particles and e^ is the single particle 



energy of the first ILL. The filhng fraction is given hj u = Np/N. Since the first ILL is 
dispersionless AEq{Np/N) is independent of t. Hence only interaction effects are taken into 
account. We search for incompressible states by computing the energy gap to adding a full 
electron, 

A,(.) = A£„ (^) - 2A£„ (f ) + AE„ {^) (6) 

as well as the energy gap to adding one flux quantum, 

A,,(.) = A£„ ( J^) - 2AE„ (§) + AE„ ( J^) . (7) 

Figure 2, plots a summary of our results for various fllling fractions and as a function of 
the lattice size. Several points may be made: i) The gap to adding a full electron shows 
more fluctuations as a function of lattice size than the gap to adding a flux quantum. 
A$„(z/ = 2/5) and A<j,q(i/ = 1/3) show a very smooth behavior as a function of lattice 
size. This provides convincing data for the occurrence of incompressible states aX v = 1/3 
and V = 2/5 on the lattice. The order of the gaps to adding a flux quantum is correct: 
A$o(i^ = 1/3) > A$o(z/ = 2/5) > A$„(z/ = 3/7) |g. Note that we obtained for particle 
excitations: Ap(z/ = 1/3) > Ap(z/ = 3/7) > Ap(z/ = 2/5) (data not presented). Finally, due 
to particle hole symmetry within the lowest ILL identical gaps at z/ = 2/3, u = 3/5, and 
u = 4/7 were observed, ii) At z/ = 1/3, and in the spirit of the Laughlin wave function, 
removing one flux quantum is equivalent to adding a quasiparticle. Since a full electron is 
composed of three quaisiparticles, one expects if the interaction between quasiparticles is 
small: A$(,(z/ = 1/3) ~ Ap(z/ = l/3)/3. This is conflrmed in flgure 2. However, it appears 
that 3A$,)(z/ = 1/3) is slightly smaller than Ap(z/ = 1/3) which support the existence of an 
attractive interaction between the quaisiparticles. Hence, the quasiparticle picture seems 
to be justifled also on the lattice, iii) The gap at z/ = 1/2 shows very large fluctuations 
as a function of system size. From our numerical data, we clearly may not obtain reliable 
information on the nature (incompressible or fermi liquid [|T^]) of the z/ = 1/2 ground state 
for the considered short range interactions. 
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Figure 3 plots the density- density correlation functions (^) for the z/ = 1/3 incompressible 
state on a 24 X 24 lattice. The lattice results (solid circles) are compared to the prediction 
of the Laughlin wave function [|T^ (solid line). For this lattice size, l^/a ~ 2 and it is 
remarkable to see that there is already nearly perfect agreement between the lattice and 
continuum results. Figure 3 also plots the quaisiparticle g{r) which is obtained by removing 
a flux quantum from the v = 1/3 state. g{r) for the quaisiparticle state is also a smooth 
function at small distances, and as expected decreases for r — ^ with a smaller power law 
than the u = 1/3 state. 

The numerical data show explicitly that the Hamiltonian Hq + Hj is scale invariant at the 
numerically accessible odd denominator filling fractions. The lattice constant is an irrelevant 
length scale, and the characteristic correlation functions scale with Iq, the only length scale 
involved. Hence, the choice of the interaction (|^) is a renormalization group fixed point. 

Mean field theories for partially filled ILL at z/ = 2m+i '^^^ t)e obtained by decomposing 
the electron operators into an odd number of fermionic partons with charges 2m+i ^ '^r ~ 
Y{a=i'4'ri subjected to the constraints '?/'°Vr = "^^r '4'r i '^Oi,(3 0. If each species fills up 
the ILL with N^ = ^m+i ^^^ mean field decoupling of the product of parton operators 
is stable if the interaction between the electrons is repulsive and short ranged ||5[. The 
parton wave functions have 2m + 1 components when we put them onto a torus. On one 
hand, this implies that the product wave function is periodic and hence is a candidate 
for an incompressible electron state 0. On the other hand the constraints on the parton 
currents induce fractional statistics [^. The important point is that in our approach the 



correct multivaluedness of anyons on a torus [iq] is automatically provided. The hierarchy 



of IQL states, which we partially observe on the lattice, can be understood analogously to 
the continuum situation. The continual condensation of localized quasiparticle excitations 
into Laughlin states relies at the mean field level on the self-similarity of the Hofstadter 
spectrum. This property assures that the levels at smaller magnetic field strengths, which 
are populated by quasiparticle excitations, have the same overall structure. The continuum 
one particle spectrum is trivially self-similar, since in this case the energies depend linearly 



on the magnetic field strength. By the introduction of the lattice the hierarchy of IQL states 
is unfold. 

To conclude, we have considered tight-binding spinless electrons on a square lattice in 
a strong magnetic field and submitted to a next-neighbor repulsive interaction. We defined 
lattice Landau levels (ILL) as well as a continuum limit which insures that they scale to 
continuum Landau levels. By projecting onto the lowest ILL and already on lattices where 
lo/a ~ 2 the FQHE at odd filling fractions was observed with remarkable precision. Those 
observations open a new region, the continuum edge of the Hofstadter spectrum, where the 
FQHE may be observed and studied. 

We are very indebted to A. Muramatsu who suggested this approach to the FQHE. We 
would like to thank W.Hanke, M. Imada, C Kiibert, K. Kusakabe, F. V. Kusmartsev and 
S. Meixner for motivating discussions. F.F.A. would like to thank the DFG for financial 
support under the grant number Ha 1537/6-2 as well as the Japanese Ministry of Education 
and S.B. the financial support of the BMFT under the grant number 03-HA3WUE. 



Figure captions 

Fig. 1 The magnetic length as a function of the lattice size. An estimate of /q on the 
lattice (symbols) is obtained by fitting the lattice value of the density-density 
correlation functions to their continuum value (^ in the case of a filled first 
landau level [0. We consider two lattice directions, r = r(l, 0) (downward 
triangles) and r = r(l, 1) (circles) as well as m = 1,3 and 5. The solid lines 



correspond to the continuum value: Iq/L = ^/m/^/2^TN where A^ corresponds to 
the linear size of the lattice. 

Fig. 2 Gaps to adding a full electron Ap(z/) (H) as well as a flux quanta A$q(z/) (0) 
for various fillings and as a function of the lattice size. The data refers to the 
special case m = 1 (see equation (0)). 

Fig. 3 Density- density correlation functions for the z/ = 1/3 state on a 24 x 24 lattice 
(solid circles). The open circles plot g{r) for the quaisiparticle on a 23 x 23 
lattice. The solid line corresponds to g{r) as obtained from the Lauglin wave 
function at z/ = 1/3 |]T^. The data refers to the special case m = 1 (see equation 
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